Precision cosmology defeats void models for acceleration 
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The suggestion that we occupy a privileged position near the centre of a large, nonlinear, and 
nearly spherical void has recently attracted much attention as an alternative to dark energy. Putting 
aside the philosophical problems with this scenario, we perform the most complete and up-to-date 
comparison with cosmological data. We use supernovae and the full cosmic microwave background 
spectrum as the basis of our analysis. We also include constraints from radial baryonic acoustic 
oscillations, the local Hubble rate, age, big bang nucleosynthesis, the Compton y-distortion, and for 
the first time include the local amplitude of matter fluctuations, ag. These all paint a consistent 
picture in which voids are in severe tension with the data. In particular, void models predict a very 
low local Hubble rate, suffer from an "old age problem", and predict much less local structure than 
is observed. 

PACS numbers: 98.80.Es, 95.36.-l-x, 98.65.Dx 
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I. INTRODUCTION 

Over the last decade the foundations of the stan- 
dard model of cosmology (SMC; see, e.g., [1]) have 
been laid. The dominant components of the SMC at 
this current epoch are baryonic and cold dark matter, 
together with dark energy, whose equation of state is 
close to that of a cosmological constant. The dark en- 
ergy is driving the accelerated expansion of a flat, ho- 
mogenous and isotropic Friedmann-Lcmaitre- Robertson- 
Walker (FLRW) background. This standard model is 
also referred to as the (flat) A cold dark matter (ACDM) 
model. 

Because of a perceived lack of elegance or naturalness 
in the SMC, several alternatives have been pursued. In 
recent years there has been a resurgence of interest in 
inhomogeneous cosmological models as an alternative to 
dark energy [2-5] (see [6] for a recent review). This has 
been partly motivated by the fact that the local expan- 
sion rate is greater in underdense regions of the Universe 
(for growing modes). If an underdensity or void is large 
enough in amplitude (density fluctuation Sp/ p ^ 1) and 
physical size (^IGpc), one can mimic the acceleration 
due to dark energy if one occupies a privileged position 
near the centre of the void. This is essentially because 
a greater expansion rate near the observer spatially can 
be hard to distinguish from an expansion rate increas- 
ing in time, for observations confined to the past light 
cone. These models break with the standard assumption 
of cosmological-scale homogeneity, and retain isotropy at 
the cost of a violation of the Copcrnican principle. 

There arc several serious physical and philosophical ob- 
jections to this scenario. First, the probability of pro- 
ducing a void of the required amplitude and size in stan- 
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dard structure formation models is extremely small [7]. 
Second, the probability of a random observer lying close 
enough to the centre of the void to avoid a large cos- 
mic microwave background (CMB) dipole is also very 
low: these models appear to strongly violate the Coper- 
nican principle [8-11]. Finally, the void would need to 
be close to spherical to match the observed isotropy in 
the Universe. These features make the model somewhat 
unappealing from the point of view of Occam's Razor. 

However, ACDM appears to suffer from severe fine- 
tuning problems of its own and it is conceivable that 
unknown physics could produce the conditions necessary 
for such a void. Therefore, it is prudent to ignore our 
philosophical prejudices and use the available cosmologi- 
cal observations to decide the issue. Given the increasing 
scope and precision of those observations, which largely 
support the SMC, demonstrating the viability of models 
that depart so dramatically from ACDM would certainly 
appear unlikely. Nevertheless, the relevance of this line of 
inquiry extends beyond particular models of dark energy, 
and refuting these inhomogeneous models would, impor- 
tantly, strengthen our confidence in the fundamental as- 
sumptions of our standard cosmological models. There 
is presently little direct evidence for cosmological-scale 
homogeneity, in particular. 

The notion that a local underdensity might be an ex- 
planation for the observed redshift-luminosity distance 
relation of Type la supernovae (SNe) appears to have 
been first proposed in Refs. [12-14]. To describe such a 
large void, the spherically symmetric Lemaitre-Tolman- 
Bondi (LTB) spacetime [15-17] is usually used. This is 
an exact solution of Einstein's equations, with a realistic 
pressureless matter source. Most early studies focused 
on fitting the SN data to LTB models of voids. This is 
not a good way of discriminating a void from ACDM, 
since it is possible to construct a radial density profile 
of a void which exactly matches the redshift-luminosity 
distance relation of ACDM (see, e.g., [18]). This high- 
lights the need to consider other observations. More re- 
cently, studies have included baryon acoustic oscillations 



(BAO) [3], spectral distortions in the CMB [2], the kine- 
matic Sunyaev-Zeldovich effect [19], and estimates of the 
Hubble rate [5]. 

Crucial evidence for the SMC is provided by measure- 
ments of the CMB, which give us a precise picture of the 
Universe at early times. The CMB also constrains the 
distance to the last scattering surface (LSS) via the an- 
gular scale of the acoustic peaks in the observed power 
spectrum. Other authors have used this angular scale by 
itself in their analysis of voids, but the radial degree of 
freedom can again be tuned to match that of ACDM [20- 
23]. 

In Zibin, Moss, and Scott [4] (hereafter ZMS), we 
stressed the importance of using the full CMB anisotropy 
power spectrum (see also [24]) and introduced the effec- 
tive Einstein-de Sitter (EdS) approach to calculating the 
anisotropics in arbitrary LTB models. This allowed us 
to place much tighter constraints on void models than 
previous studies. We also explored a very wide region 
of void radial profile space using a spline parameteriza- 
tion, which allowed us to identify models which fitted 
the CMB + SN data. However, we found that models 
which matched the CMB power spectrum had a local 
Hubble rate. Ho, so low as to rule them out. We ar- 
gued that it would require significant fine tuning of the 
primordial perturbation spectrum in order to circumvent 
this conclusion, but we did not investigate this in detail. 
More recently, Ref. [25] stressed the importance of ob- 
taining model-independent constraints from the CMB in 
a broader context. 

Another class of observations that could potentially 
put strong constraints on void models of acceleration is 
that of structure in the late Universe. The evolution of 
matter perturbations on LTB backgrounds is expected to 
differ significantly from that in the SMC, so observations 
of the late-time matter power spectrum shape or am- 
plitude, weak lensing, or integrated Sachs- Wolfe (ISW) 
effect would likely help to rule out void models. Unfor- 
tunately, the evolution of perturbations on LTB space- 
times is an extremely difficult problem [26] , and no exact 
calculations have yet been carried out. An approxima- 
tion method which drastically simplifies the calculations 
was proposed in Ref. [27] (see also [28]); however, again, 
no confrontation with data has yet been performed. A 
large part of the difficulty is that observations, e.g. of the 
matter power spectrum, are often presented in a model- 
dependent manner, in that an underlying FLRW back- 
ground is assumed. 

In this work we undertake an extended update of ZMS, 
using the latest (and additional) cosmological data. We 
improve our method in several respects, and also con- 
sider extensions to the "basic" void model. We begin 
in Sec. II with a summary of the exact general relativis- 
tic framework within which we perform our analysis. In 
Sec. HI we describe our effective EdS approach to cal- 
culating CMB spectra, and describe the other data we 
consider, as well as our statistical approach. We also 
enumerate our improvements to the approach of ZMS, 



which include updated CMB and supernova data and 
advances in our numerical and fitting procedures. Our 
analysis also includes new types of observations, such as 
the y-distortion of the CMB frequency spectrum, as well 
as the first constraints on void models using the ampli- 
tude of local structure. We avoid the problem of model 
dependence in examining structure data by considering 
only the most local observations. Section III closes with 
the constraints which form our main results. We then 
consider whether modifications to the primordial spec- 
trum can alleviate the problems facing voids in Sec. IV, 
before examining an interesting multi-valued region of 
parameter space. We close with conclusions in Sec. V. 
Throughout this paper we set c = 1. 



II. OVERVIEW OF LTB FRAMEWORK 

A. Exact solution 

We consider an observer in an isotropic, inhomoge- 
neous Universe which at late times is dominated by pres- 
sureless matter. Such a model is described by the LTB 
spacetime, in which the Einstein equations can be solved 
exactly. The line element of the LTB metric is given by 



ds^ = -rfr 



y- 



1- A' 



-dr' + Y^Q' 



(1) 



where a prime denotes the derivative with respect to 
comoving radial coordinate r, and t is the proper time 
along the comoving worldlines. The curvature function 
K = K{r) is a free radial function (bounded by K < 1), 
and the areal radius Y = Y{t, r) is given parametrically 
by the exact solution 
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for K — 0. Here, the free radial function te = ^b(?') is 
known as the "bang time", since the cosmological sin- 
gularity occurs at t = tBir). In this work we assume a 
homogeneous big bang and set <b = 0. In this case the 
void only has a growing mode [27, 29], and hence we can 
use standard cosmological perturbation theory up to the 
epoch where the void becomes nonlinear. The function 
M = M(r) is a further arbitrary radial function. How- 
ever, one of these free radial functions is a gauge mode 
corresponding to a rescaling of the radial coordinate. We 
set M = r^ , which imposes no restriction on the solutions 



except that they wiU not be vahd past the "equator" of 
a spatially closed model. The LTB spacetime is then 
completely specified once the single free function K{r) is 
fixed. 

Covariant physical quantities within the LTB space- 
time are given by (see, e.g., [27]) 
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Hn + 2H-Y , 
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where p, 6, and E are the comoving matter density, ex- 
pansion, and shear, respectively, and ^'^^i? is the Ricci cur- 
vature of the spatial comoving-orthogonal hypersurfaces. 
The radial and transverse expansion rates are given by 
Hr = Y'/Y' and Ht = Y/Y, where the ovcrdot de- 
notes the derivative with respect to t. We define the local 
density parameter by fij^'^ = 247rGp/0^. These quanti- 
ties reduce to the standard FLRW expressions {9 = 3H, 
S = 0, fi|,°'= = Qm) when Y{t, r) = a{t)r, where a{t) is the 
FLRW scale factor and H = a/ a the Hubble rate. Along 
the centre of symmetry in an arbitrary LTB spacetime, 
we have S = and TJr = iJx = H . 

The null radial incoming geodesies [ds^ = dVL^ = 0) 
are described by the equations 



dt 
dr 



-Y' 



dz 
dr 



{l + z)Y' 



(9) 



where the redshift z is measured along the past light cone, 
increasing from z = at the origin. The luminosity and 
angular diameter distances from the centre to redshift z 
are then 



dL = {l + zyY, dA = Y. 



(10) 



The LTB equations were coded in a Fortran module 
to interface with the CMB code CAMB [30]. We first 
relate the curvature function to the early density profile 
by Eq. (11), and discard models with 1 — K{r) < 10"'^ to 
ensure numerical stability. Given the Hubble rate Hq = 
100/iokms~^ Mpc"^ observed today (t = to) at the void 
centre, we evaluate the time to using the exact solutions. 

We then integrate Eq. (9) along the past light cone 
from the centre today to a redshift far outside the void. 
During integration, we store quantities such as t, z, and 
c^L with the integration variable r in a finely-spaced ar- 
ray. These are then interpolated using cubic splines when 
intermediate values are required (such as fitting redshift- 
luminosity distances of SNe). 

In some models (typically very deep voids) the ra- 
dial expansion rate can become negative, when overdense 
regions break from the background expansion and be- 
gin to contract. As can be seen from Eq. (9), negative 
Hn = Y'/Y' implies that the redshift will decrease with 
increasing radial coordinate down the light cone. This 
leads to an interesting class of models with multi-valued 
distance-redshift relations [31] (see also [32]). In extreme 
cases the redshift can even become negative. We flag 
these multi-valued models and do not use them in our 
basic constraints, due to uncertainties in interpreting fits 
to SNe. The multi-valued cases will be discussed further 
in Sec. IV B. 

After overdense regions contract sufficiently, they can 
experience a shell- crossing singularity, when Y' crosses 
zero and hence the density diverges according to Eq. (5). 
These singularities are a symptom of the pressureless as- 
sumption of LTB: in reality, interactions on small scales 
would modify the dynamics and the LTB solution would 
become invalid. Therefore, we perform checks in our code 
and discard models in which Y' = anywhere on the past 
light cone. (Of course, shell crossings may still occur in 
the future, but we do not need the solution in such un- 
obscrvablc regions.) 



B. Numerical implementation 

In our numerical calculations, instead of fixing K{r) we 
found it more convenient to express the single free radial 
function by the early-time comoving density perturba- 
tion, 6{ti, r) = [p(ti, r) - pFLRw(^i)] /PFLRw(ij): where ti 
is a time near last scattering. Because we only consider 
the growing mode of the void, S{ti,r) can be considered 
a linear fluctuation from FLRW. In this case the curva- 
ture function K{r) can be determined from the density 
perturbation using the linear growing mode relation be- 
tween the comoving density perturbation and curvature 
(see, e.g., [27]) and Eq. (8). The resuh is 



{rxy 



AOttG 



PF-LRw{ti)Y^S{ti,r) , 



fll) 



where in the early linear regime Y — a.^r and a^ ~ a,(ti) 
is the initial scale factor, which we arbitrarily set to a.; = 
10-3. 



III. BASIC CONSTRAINTS 

In this section we detail the basic cosmological con- 
straints on LTB models. We have made several changes 
and additions to our analysis in ZMS, most of which are 
discussed in more detail later. In summary, these arc 

• Cosmological data have moved on. We use CMB 
data from the Wilkinson Microwave Anisotropy 
Probe (WMAP) 7-ycar release [33, 34] and the lat- 
est Union2 SN compilation [35], which has nearly 
double the number of SNe and improved outlier re- 
jection compared with our previous analysis. 

• We change the integration variable of the LTB 
equations to use the radial coordinate rather than 
redshift. This allows us to find models with multi- 
valued distance-redshift relations and to more eas- 
ily reject models with shell crossings. 



• A more efficient metfiod is used to searcfi tfirougli 
the void parameter space wlien fitting to CMB ob- 
servations. This involves sampling from a set of 
effective parameters derived from CMB data. 

• We investigate an alternative parameterization of 
the void profile, which allows us to separate the 
physics of fitting voids to SN and CMB data. 

• We investigate voids embedded in a spatially 
curved FLRW background. 

• We discuss additional cosmological constraints, 
such as those from the Compton j/-distortion of 
the CMB blackbody spectrum. We also present 
the first estimates of the local amplitude of matter 
fluctuations smoothed over 8/iQ^Mpc spheres, cts- 



A. Formalism 

1. CMB spectra 

The effective EdS method introduced in ZMS forms the 
basis for our computations of CMB power spectra in LTB 
models. The method can also be applied to any model 
with an expansion history different from EdS (FLRW 
with pure dust source and vanishing spatial curvature) 
from recombination to today. The basic idea is as fol- 
lows: we find the parameters of an effective EdS model 
which has the same physics at recombination as the void 
model, as well as the same angular diameter distance to 
the LSS, in proper units at last scattering. This ensures 
that the Q's will be identical between the effective EdS 
and LTB models, apart from any sources of secondary 
anisotropy between today and the LSS. The spectra can 
then be readily calculated by feeding the effective param- 
eters into public CMB codes. Using EdS for the effective 
model is convenient since it ensures that no secondary 
ISW component will be present. 

To generate the effective model, we must match the 
angular diameter distance to the LSS between the EdS 
model and the void model, and also match the Hubble 
rate and density of relativistic species and matter com- 
ponents at the LSS (i.e. the physics at last scattering). 
The extra degree of freedom required in the effective EdS 
model comes from our ability to specify a CMB temper- 
ature for the effective model, different from the actual 
To = 2.726 ± 0.001 K [36] which we observe today. (Note 
that the effective model will include substantial radia- 
tion at early times, and hence should not properly be 
called "EdS" . However, for the purposes of feeding effec- 
tive late-time parameters into CMB codes, it is essen- 
tially EdS.) 

We follow the procedure given in ZMS, generalized 
slightly to accomodate void models that asymptote to 
spatially curved FLRW, rather than just to flat FLRW. 
To specify an LTB model we must fix the free radial pro- 
file and the Hubble rate at the observation point today. 



Hq. Once this is done, we first compute the coordinates 
(^1117 ''m) a-t a "midpoint" redshift Zm down the light cone, 
far outside the void. This redshift is chosen such that 
the background LTB shear is negligible, and the radi- 
ation density (treating radiation as a test field on the 
LTB background) is small compared to the matter den- 
sity. In our numerical calculations we use Zm = 100, 
and have checked that our results arc not affected sig- 
nificantly within the range 50 < Zm < 200. We cannot 
choose Zin to coincide with the LSS, since radiation is 
important at background level there and the LTB solu- 
tion does not include radiation. However, as long as the 
background shear is negligible beyond Zm, matching at 
Zm will be essentially equivalent to matching at the LSS. 
The basic idea is that where radiation is unimportant (for 
z < Zm), the LTB model can be used, while where back- 
ground shear is unimportant (for Zi„ < z < zlss)- the 
true inhomogeneous spacetime (i.e. including radiation 
as a source) is very closely approximated by an FLRW 
matter plus radiation model, so that matching at z,n is 
essentially equivalent to matching at zlss- 

We then compute the Hubble rate, i/ni, at Zm in the 
LTB model and integrate back up the light cone into the 
effective EdS model to comoving coordinate r^*^^ = 0. In 
the case that the void becomes asymptotically spatially 
flat, we ensure that the proper distance to the LSS is 
correct by setting a{r^^^)r^^ = Y{tn^, r^), for EdS scale 
factor a. This allows us to calculate the effective EdS 
mean temperature and Hubble rate via 



\ i -r Zi^ 



where 



1 + Z^dS 



l + z„ 



(1 + aEdS^Eds^^^/2)^ 



(12) 
(13) 



(Note that the quantity z^'^^ was called z^^ in ZMS, and 
that in general z^'^^ ^ 0, since we choose the values of 
the rcdshifts Zm at the "midpoint" to be identical in both 
void and effective EdS models.) The effective parameters 
Tff '^^ and H^'^^ define the point in EdS which observes 
the same primary CMB as the specified LTB model. 

In the case that the LTB model asymptotes to spa- 
tially curved FLRW, we still choose our effective model 
to be EdS, but we must modify the above calcula- 
tion. In particular, we modify the matching condition 
to a(rEdS),.Eds ^ Y{t^,r^) + Ac^a, where 



Ada 



Ma(^lss) 
Ma(^lss) 



^A l^m) I curved 
FLRW 

dA{Zm)\ flat 
FLRW 



(14) 



is the difference in the angular diameter distance incre- 
ment from Zni to Zlss between the actual curved FLRW 
model and flat FLRW. In this expression, the curved 
FLRW values are calculated in the FLRW model which 
has the same matter and radiation densities and curva- 
ture at z,n as does the actual void model. For flat FLRW, 



the increment from Zm to zlss is calculated directly by 
integrating the relation 



dr 



dz (1 



1 



.)Hizy 



(15) 



For the curved FLRW case, Eq. (15) only provides the 
coordinate increment amAr; this must be further trans- 
lated into an increment in c^a using the relations 



1 



i/„V|l-f^m^(^m)r 



dA 



(16) 



(and similarly for the open FLRW case). With this cor- 
rection, the effective EdS model will produce essentially 
exactly the correct primary CMB as long as the LTB 
model becomes essentially homogeneous (shear-free) by 



It is easy to evaluate the effective EdS parameters for 
a ACDM cosmology instead of an LTB model. They can 
be found by substituting the quantities 



Hm ~ -ffoV ^m (1 



,EdS 




into Eq. (12), where fl^ is the matter density param- 
eter today in ACDM. The integral in Eq. (17) can be 
computed numerically, and is insensitive to the choice 
of Zjn at the level of < 0.1% for z,„ > 5. Plots of the 
effective parameters are shown in Fig. 1. As 0,n — > 1 
we find T^'^^ -^ Tq and H^'^^ -^ Hq as expected, while 
for low Onj the effective parameters are very different 
from the actual observed parameters. For the WMAP 
7-year best-fit ACDM model, given by J7m = 0.262 and 
Ho = 71.4kms-iMpc-i [34], we find Tffds ^ 3.416K 
and H^'^^ = 51.2kms~^ Mpc~^. A void model must 
have effective EdS parameters close to these values if it is 
to produce CMB power spectra similar to those actually 
observed. 



2. Low I CMB 

The Cf's of the effective model are essentially identi- 
cal to the actual model on intermediate and small angu- 
lar scales. For temperature, however, the ACDM model 
creates an ISW anisotropy on large scales, due to decay- 
ing gravitational potentials along the line of sight. A 
void model is also expected to create an ISW signal, but 
the complications of evolving perturbations on an LTB 
background have prohibited a rigourous calculation of the 
LTB ISW effect, so that the size and shape of the ISW 
signal is uncertain. 

In order for the effective EdS method to remain viable, 
however, we introduce a marginalization uncertainty over 
an ISW template; large angle data are necessary for con- 
straining the amplitude and spectral index of initial fiuc- 
tuations, so one would lose important predictive power 
by applying a simple cutoff in £. Specifically, the Ci spec- 
trum is given by 



where A£(fc,io) is the associated multipole moment for 
the photon distribution and Vsik) is the initial comoving 
curvature power spectrum, parameterized by Vsik) = 
Aq, (k/ko) ^ , where As is the initial scalar amplitude, 
TT-s the scalar spectral index, and fco the pivot scale (which 
we fix to 0.05 Mpc~^). We modify the transfer function 
to 



Af (fc, to) = A^^^(fc, to) + Aisw^'r (k, to) 



(19) 



a = 4TT d{hik)Vsik)\A,{k,to)f 



(18) 



where Aisw is the ISW amplitude. Since the data on 
large scales are limited by cosmic variance, the precise 
shape of Ap^(A;, tq) is not very well constrained, and so 
we use a template from the WMAP7 ACDM model. This 
choice provides a conservative estimate of the expected 
ISW signal from a void. 

The effective EdS parameters can now be used in any of 
various CMB anisotropy codes; we used CAMB [30]. The 
full set of parameters specifying the effective model are 
the baryonic matter fraction /b = pb/pm (which we take 
to be spatially constant throughout this work); Tif"^^; 
i/Eds ^ ioo/i^dS]^ms-iMpc-i; As; ng; Aisw; and the 
redshift of reionization zf^^. We show the temperature 
and polarization power spectra of the effective EdS model 
for WMAP7 ACDM in Fig. 2. Here, we have fixed the 
reionization optical depth to the WMAP7 value, so that 
Zro is lower in EdS compared to ACDM (~8 as opposed 
to -10). 

For the TT spectrum, the agreement between Cg's is 
excellent at all i. For polarization, particularly the TE 
cross spectrum, there is still a discrepancy on large scales. 
These modes are sourced at z < Zrc, so presumably the 
difference arises from the slightly different reionization 
history of EdS and (for TE) the cross-correlation with 
our simplified ISW template. 

The generation of large angle polarization is an inter- 
esting question in LTB models. Electron scatterers in 
the region of the void at z < z^-c will see a tempera- 
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FIG. 1: Effective EdS parameters T^'^^ and H^"^^ as a func- 
tion of the matter density Qm in the corresponding ACDM 
cosmology. 



ture anisotropy dominated by a dipole (we discuss this 
further in the following section). Each scatterer will also 
see a smaller local quadrupole, with a relative value larger 
in the void periphery. This quadrupole does not induce 
any polarization for an observer at the centre of the void 
though, since it is parallel to the radial direction. 

The void will, however, have a different reionization 
history than for the effective EdS model. Our earlier 
work in ZMS showed that a large overdense outer shell 
at z > I is required to fit the CMB data (in the absence of 
background spatial curvature), and the precise details of 
this shell will modify the reionization history, and hence 
the low £ polarization. For this reason, together with 
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FIG. 2: CMB power spectra for the WMAP7 ACDM cos- 
mology (solid curves) along with the corresponding effective 
EdS model (dotted curves). We show the TT (top panel), EE 
(middle panel), and TE (bottom panel) spectra. 



our simplified treatment of the ISW cross-correlation, we 
choose to ignore TE data on large scales, using a cutoff of 
^ = 50. Since Zrc is less constrained by the data without 
large scale TE, we apply a conservative prior on z^^ in 
our analysis, the details of which are given in Sec. fllB. 



3. Dipole anisotropy 

Off-centre observers in an LTB model will see a CMB 
dipole, AT/T = (3 cos 9, where d is the angle with respect 
to the radial direction. This dipole will be the dominant 
anisotropy sufficiently close to the centre. As we will see, 
even for central observers this will have important con- 
sequences: the scatterers along the observer's light cone 
will produce spectral distortions in the CMB seen at the 
centre. To compute f3{z), the dipole seen by scatterers 
at redshift z down the light cone of the central observer, 
we can consider the propagation of incoming and outgo- 
ing radial geodesies towards the scattcrer, and assume 
the anisotropy is a pure dipole. (Incoming/outgoing re- 
fer to the direction at the scatterer.) More detailed cal- 
culations using non-radial geodesies show that this is a 
good approximation, apart from in the peripheral region 
of the void and further outwards, where the anisotropy 
subtends a smaller angle on the sky [8, 37]. 

To calculate /3(z), we follow a similar procedure to cal- 
culating the effective EdS CMB parameters. Concretely, 
we first evaluate the coordinates (tg, Vg) of a scatterer 
at Zg. We then continue the light cone integration to a 
position far from the void centre, choosing Zin = 1000 
(i.e. approximately the redshift of the LSS), to find the 
originating coordinates (iin,nn) of the radial incoming 
geodesic. To complete the calculation, we compute the 
redshift Zout of an outgoing geodesic, also originating at 
time tout = tin, which propagates back through the void 
centre and arrives at the same {ts,rs) as the incoming 



J 



geodesic. We use a value of Zin = 1000, larger than 
Zni = 100 used in the CMB calculations, since the outgo- 
ing ray generally passes through the void and the spatial 
curvature (and shear) is non-neglible for some models 
within the void at t{z = 100). Moreover, the precise ori- 
gin of the photons can affect the dipole at the observer, 
as shown in Sec. Ill C 3, so we choose the LSS explicitly. 
The condition tout = tin ensures that both radial rays 
originated at essentially the same density, which defines 
the moment of last scattering in the absence of isocur- 
vature modes. The dipole anisotropy observed at Zg is 
then 



^out Z\yi 



Z, 



(20) 



out ~r ^in 



This dipole is indistinguishable from a peculiar velocity 
at the scatterer. 



4- Compton y-distortion 

Non-local processes can also provide complimentary 
constraints on LTB models [2, 38]. These arise from the 
re-scattering of photons in the region of the void after 
the Universe reionized at z^^. One such process is the 
Compton y-distortion of the CMB blackbody spectrum. 
The j/-distortion in the single-scattering and linear ap- 
proximations is given by [2, 39] 



2/(n) 



3 

32^ 



dz 



dr 
dz 



d^h' 



l + (n-fi')^ [e(n',n',2) 



e(n,n'. 



(21) 



r 



where 8 = AT/T is the temperature anisotropy at the 
scatterer and r is the optical depth. The vector fi is the 
direction from observer to scatterer and n' the direction 
of the photon arriving at the scattcrer from the LSS. 
When the temperature anisotropy is dominated by the 
induced dipole, we have 0(n, n',z) — (5 (z) cos 9, with 
cos9 = n • fi'. 

When the dipole dominates, we can complete the an- 
gular part of Eq. (21) analytically to leave 
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dr 



y = T7^ dz-rP{zr , 



dz 



(22) 



where we integrate up to the redshift of reionization z^c . 
The redshift dependence of the optical depth is given by 



dr 
dz 



aTTieiz) 



dt 

dz 



cJT9^h{l~Yno/1)^t'{z)dt 



247rGm„ 



dz' 



(23) 



with (Tt being the Thomson cross section, rig the electron 
number density, Ypje the helium mass fraction, and nip 
the proton mass. In our calculations we fix Yhc = 0.24. 



5. Baryon acoustic oscillation scale 

LTB models generically exhibit shear at background 
level, so that H^ ^ Ht for r > 0, and this anisotropy in 
the expansion increases as the void grows and becomes 
nonlinear. Therefore, proper length scales, which are 
isotropic in the FLRW regime at early times, also be- 
come different in the radial and transverse directions at 
late times. One scale of particular cosmological inter- 
est is the sound horizon at the drag epoch, defined as 
the time when baryons are released from the Compton 
drag of photons. This can be approximated in an FLRW 



cosmology by [40] 
Spiz) 



44.51n[9.83/(an^o) 



{i + z)^i + To{nJ^' 



:MpC. 



(24) 



Importantly, Sp{z) is the proper sound horizon at red- 
shift z, which is a model-independent measure. In Fourier 
space, this scale leads to the characteristic series of 
baryon acoustic oscillations observed in the matter power 
spectrum. 

To calculate the observed BAO scales in the LTB 
model, we first generate a new effective EdS model, fol- 
lowing the method introduced in ZMS. This effective 
model will share the same early physics as the LTB 
model, but unlike the case of the CMB, here there is no 
need to match the angular diameter distance to the LSS. 
In order to allow us to use Eq. (24), we choose T^^^ ~ Tq 
for the effective model. Substituting z^'^^ = in Eq. (12) 
then gives 



H, 



EdS 



H„ 



(1 + 2m)3/2 ' 



(25) 



where we again choose Zm = 100. To evaluate Sp(zm), we 
then substitute this effective Hubble rate, together with 
r^in = 1 and rib = /b, into Eq. (24), checking that these 
effective parameters are within the regime of validity of 
that approximation. 

Since the effective EdS and LTB models share the 
same early physics, we can conclude that the (essentially 
isotropic) sound horizon at coordinates (f(zm),r(zm)) = 
(tm, rm) in the LTB model is given by Sp(zm). But, at the 
early time tm, the LTB spacetime is very close to FLRW. 
Therefore, we can conclude that the sound horizon takes 
essentially the same value Sp{zm) at (t^^r), for any r 
in the LTB model (assuming no significant isocurvature 
mode associated with the void). Thus, to evaluate the 
BAO scales at redshift z down the observer's light cone, 
i.e. at coordinates (i(z), r{z))^ we only need to evolve the 
scale Sp(zm) according to the radial and transverse linear 
expansions in the sheared LTB background between the 
points (t„i,r(z)) and (t(z),r(z)). The radial and trans- 
verse proper BAO scales at redshift z are then, respec- 
tively. 
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(26) 
(27) 



Finally, we can rewrite the radial and transverse proper 
BAO length scales as corresponding redshift and angular 
intervals at redshift z in the LTB model. The result is 

Az{z) - il + z)l^^°iz)Hnit{zlriz)), (28) 



Ae'(z) 



/BAO(^) 

Y{t{z),r{z))- 



(29) 



As stressed in ZMS, since these final quantities are 
directly measurable, they are model-independent, and 



hence can be unambiguously compared with data. In ad- 
dition, the radial scale Az{z) is expected to be a strong 
discriminator of void models. This is because in the 
peripheral region of the void the radial expansion rate 
HYi{z) is generically suppressed, as dust from the void 
piles up. This in turn leads to a suppression in the ra- 
dial BAO length scale /{^^^(z) in that region. Therefore, 
according to Eq. (28), these two suppressions reinforce 
each other, resulting in a heavily suppressed radial red- 
shift interval Az{z) in the void periphery. 

In the literature on BAO observations, constraints are 
often expressed for the isotropized distance measure (see, 
e.g., [41]) 



Dviz) = 



.W 



Hiz) 



1/3 



(30) 



We do not attempt to estimate Dv{z) for the LTB mod- 
els for two reasons. First, comparison with Eqs. (28) and 
(29) reveals that Dy is related to a combination of radial 
and angular scales, but weighted more heavily to the an- 
gular scale. The angular scale is determined in part by 
the angular diameter distance, d^iz) = Y(z), which, for 
an LTB model that fits the SN data, must be close to that 
of the standard A model. Therefore we do not expect Dy 
to be a strong discriminator of void models. The second 
reason we avoid Dy is that in its construction it assumes 
an FLRW background, and hence constraints on Dy are 
not model-independent. 



6. Amplitude of local matter fluctuations, as 

The evolution of structure on general LTB back- 
grounds is a very complex problem (see [26], and [27] 
for an approximate approach). Therefore, no predictions 
for the matter power spectrum in LTB models has yet 
been made. This is unfortunate, since considerable in- 
formation is available to us in the form of measurements 
of the power spectrum, and there are likely to be signifi- 
cant differences in the evolution of perturbations between 
LTB and FLRW models. 

However, in a certain region of the LTB spacetime we 
can in fact accurately calculate the evolution of perturba- 
tions, knowing only how the structures evolve on FLRW 
backgrounds. In Ref. [27] it was shown that the comov- 
ing matter density and expansion perturbations evolve 
according to 



6p 



iirGSp - 3E(5S, 



(31) 
(32) 



where STi is the shear perturbation associated with the 
structure. These equations apply in general regions of the 
LTB spacetime. However, they are not closed because of 
the (51] term, and that term will generally couple to per- 
turbations in vector and tensor modes, leading to the 
immense complexity of the general problem. But at the 



origin we have S = by symmetry, and so the scalars 5p 
and 50 decouple from vectors and tensors. Similarly, close 
enough to the origin that E/6' ^ 1, we expect the decou- 
pling to persist to good approximation. Note that even 
though tensors generated where E is large may propagate 
to the origin, near there they will not affect the evolution 
of the scalars at linear order because of the decoupling. 

Therefore, where we can ignore the shear coupling term 
to a good approximation, Eqs. (31) and (32) close, and 
their solutions are identical to those in the FLRW case. 
Similarly, the background CYohiiion (pit), 9{t)) of the LTB 
spacetime near the origin will match that of a spatially 
curved FLRW model. The evolution of the scalar pertur- 
bations near the origin will then be the same as that in 
the FLRW model with the same background evolution as 
the origin of the LTB model. It is simple to verify that 
for an LTB model with local Hubble rate Hq and density 
parameter il\°^ at proper time to at the origin, the cor- 
responding FLRW model would have the same Hubble 
rate and density parameter 0^ = i^J^'^ (and vanishing A) 
at the same tg. For the (undcrdcnsc) void models con- 
sidered here, this corresponds to an open FLRW model. 

To calculate the matter power spectrum, we again 
use an effective model approach. As before, the effec- 
tive model will share the same early physics as the LTB 
model, but in this case the effective model will be the 
open FLRW model with parameters Hq^°^^ — Hq and 
^opcn ^ ^ioc(^ ^0). To determine the effective temper- 
ature, we first integrate down the LTB observer's light 
cone to redshift z^ to find t^ = t{zm)- Then we evalu- 
ate the Hubble rate He and density parameter ^^^% at 
coordinates (t,„, 0), i.e. at the centre at im, using the ex- 
act LTB solution. Then the temperature in the effective 
model can be calculated using the FLRW relation 



B. Method 
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(33) 



To calculate the matter power spectrum (or its ampli- 
tude on scale 8//io Mpc, i.e. as) near the origin in the 
LTB model, we then feed the effective parameters Hq^"^, 
f^°P™, and ro°P™ to CAMB. Of course the result wiU de- 
pend on the primordial amplitude of perturbations, Ag, 
which is an input parameter to CAMB, but this amplitude 
will be constrained when we fit to the CMB data. Since 
we are assuming pure growing mode LTB models, the 
same primordial amplitude should apply to early times 
at the centre as at the LSS. 

Note that in general wc will have Tq^'^'^ ^ Tq, i.e. the 
effective model will have a different radiation density to- 
day than the LTB model. However, the radiation density 
is so low today that it has negligible effect on the evolu- 
tion of matter perturbations. At early times, when radi- 
ation is important, our effective model shares the same 
physics as the void model, so the evolution of matter 
perturbations will be the same in the two models at all 
times. 



To calculate the CMB spectra for each void model we 
only need to know the effective parameters Tjf'*^, Hq'^^ 
for the model (together with the primordial spectrum 
parameters, of course). Since these effective parameters 
can be calculated very quickly, while the corresponding 
CMB spectra take considerably more time to generate, 
it is worthwhile to first investigate the posterior likeli- 
hood of these parameters (along with the other "stan- 
dard" cosmological parameters) when applied to CMB 
data. If the posterior can accurately be characterized 
by a multi-variate Gaussian around the maximum likeli- 
hood, efficient searches of the void parameter space can 
be achieved by sampling from the covariance matrix of 
the effective parameters, rather than by re-fitting CMB 
data for each void model. 

We use COSMOMC [42] to generate Markov-Chain- 
Monte-Carlo (MCMC) chains to estimate confidence lim- 
its on the effective parameters. We use CMB data 
from WMAP7 [33, 34], together with those from the 
ACBAR [43], Boomerang [44], CBI [45], and QUaD [46] 
experiments, which observe to higher i. To compare 
the goodness-of-fit, we first fit a "vanilla" ACDM model 
to our modified likelihood routine (ignoring polarization 
data at ^ < 50). The parameters in this model are: the 
baryon density, f^b^o! '^^^'^ dark matter density, flc^oi ^o; 
Zro; ^s; and ng. We marginalize over the SZ amplitude 
Asz, assuming the Komatsu and Seljak template [47], 
and apply a prior of Zj-o > 8. We do not include lensing 
of the CMB in our analysis, since at the moment we do 
not fully understand how LTB models modify the lensing 
signal. 

For the effective EdS model, one degree of freedom 
is removed since J7m = 1, but we fit for Tq'^^ and 
marginalize over Ajsw- We apply a conservative prior 
of 5 < zf^^ < 15, since the void model could have a 
lower effective z^^ than the actual z^^. In addition, we 
also investigate the case of spectral index running, such 
that the initial power spectrum is characterized by 

ln7's(fc) = ln^s + (ns - l)ln (fc/fco) + ^ [ln(fc/fco)]' ■ 

(34) 
In Fig. 3 we show a selection of 2D likelihoods from 
the MCMC chains for the effective EdS model param- 
eter space. We define the best-fit likelihood relative to 
ACDM as Ax^ = — 21og(£A/£Eds), such that positive 
values favour EdS. We find Ax^ = 0.6 for the case of no 
running and 2.5 with running. The marginalized values of 
TpEdS ^ (3 43 ^ 0.08) K and h^'^^ = 0.512 ± 0.008 (with- 
out running) are consistent with the best-fit WMAP7 
ACDM effective parameters calculated in Sec. IIIAl. 
We note that even with running, T^'^^ > 3.17 K at 2cr. 
This is important for void models, since an overdense 
outer shell (or positive spatial curvature) is required for 
yEds ^ To, as highlighted in ZMS. Without such a fea- 
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ture, one would have to consider non-trivial modifications 
of the primordial spectrum to fit the CMB if T^'^^ w Tq. 
We return to this issue in Sec. IV A. 



2. Void profile 

We define the initial radial density profile S{ti,r) of 
the LTB model in three ways. In the first, we repeat the 
analysis carried out in ZMS with updated cosmological 
data. Here, we fit a three-point cubic spline to the initial 
density fluctuation 6j = S{ti,rj), where j = 1, 2, 3. We 
fix ri = and enforce the void to be smooth at the 
origin and to smoothly approach EdS at large r by setting 
S'{ri) = (^'(ra) = S3 = 0. There are thus a total of four 
free profile parameters, and we find that additional spline 



points do not significantly improve the fit to the data. 

In ZMS it was noted that a large overdense shell out- 
side of the central void is required to fit the CMB power 
spectrum in asymptotically flat models. For the three- 
point spline parameterization, this shell will smoothly 
transition from the underdensity near the void centre. 
For this reason, we also find it convenient to decompose 
the initial density perturbation into the linear combina- 
tion of profiles 



6iU,r)^Suir)+Soir). 



(35) 



where Su represents the initial central underdensity and 
60 the initial overdense shell. These profiles are given by 
the polynomial functions 
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(37) 



such that S{ti,r) smoothly matches to an EdS region be- can be obtained by setting 

tween the void and shell. There are a total of five free 

parameters: the width and depth of the void, rg and Sq; r _ 

and the width, height, and position of the shell, Ar, 5i, 

and ri. However, the fit to the CMB power spectrum. 



~Sor^ 



Ar(15r2 + 2Ar2) 



(38) 



requiring Tq'^^ « 3.4 K, is primarily a function of the 
integrated shell density, so in our analysis we arbitrarily 
fix the shell position and width and vary only the height. 
We choose the position so that the shell lies beyond the 
SN data (z > 1.5). 

This "polynomial + shell" decomposition is useful as 
it separates the two main physical effects: the fit of the 
SN data to the central void; and the fit of the CMB to 
the outer shell. We find that some additional constraints 
(such as the Compton y-distortion and radial BAO scale) 
are sensitive to the details of the shell, which can more 
easily be investigated without the constraints the spline 
imposes on the shell position. 

In ZMS, we made the distinction between "con- 
strained" profiles, for which V = J Sp{ti,r)r^dr < 0, 
and "unconstrained" , which are free. For the polynomial 
-|- shell model, a perfectly compensated void with 1^ = 



We do not enforce this condition. Therefore, although 
our early-time profiles always satisfy Sp(ti,r) = for r 
lying between void and shell or exceeding some maximum 
radius, for the general uncompensated case the spacetime 
will depart from EdS at late times in these regions, due 
to the gravitational effect of the inhomogeneities. 

Finally, we also consider voids embedded in a spatially 
curved FLRW background. Here, the curvature required 
to fit the CMB is distributed homogeneously rather than 
confined to a large overdense shell. In this case we set 



S{U,r) =Su{r)+Si 



(39) 



where the central underdensity Su{r) is given by the 
polynomial function (36) and Si is a constant. For this 
"polynomial -f curvature" parameterization, the curva- 
ture function at large r then takes the FLRW form 
K{r) ex r^. This corresponds to a spatially constant 
density parameter fl\°'^ in this asymptotic region. Note, 
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FIG. 3: Constraints on the effective EdS parameters from tfie CMB power spectrum data. Likelifiood contours show the 68% 
and 95% confidence levels. Dotted and solid contours are for models with and without spectral index running, respectively. 



however, that there is ambiguity in choosing the time 
at which to specify such a parameter in a globally non- 
FLRW background, and observations do not directly con- 
strain the asymptotic curvature. 



3. Data fitting 

Measurements of the apparent magnitudes of Type la 
supernovae are important for constraining the void pro- 
file. In this work we use the recent Union2 compila- 
tion [35] of 557 SNe in the range z = 0.015-1.4, which is 
nearly double the number of SNe used in ZMS, with im- 
proved rejection of outliers. The residual fi between the 
apparent and absolute magnitudes, m and M, is defined 
as 



M 



Af = 51og 



10 



Mpc 



25. 



(40) 



In our fitting we adopt the standard procedure of 
marginalizing over the unknown absolute magnitude. 

For the CMB we tested the spline profile both by sam- 
pling from the covariance matrix of the effective param- 
eters and by re-fitting the CMB data each time. The re- 
sults in both cases were similar, but the former method 
was much faster, with combined CMB + SN constraints 
taking only a few hours to obtain. In our subsequent 
analysis we therefore use the covariance sampling. 

As in ZMS, we applied a conservative prior of il^°^'^ > 
0.1 at the void centre. This is consistent with estimates 
of the minimum matter density in the local Universe [48, 



49]. We investigate models with lower ^^^'^(z = 0) 
Sec. IV B. 



C. Results 

1. Basic parameters 

In Fig. 4 we plot the local density parameter il\°'^{z) 
along the observer's past light cone for several spline, 
polynomial -f- shell, and polynomial -f- curvature profiles. 
These are sampled from the MCMC chains, using CMB 
+ SN data, with the grayscale level indicating the relative 
likelihood. Each curve has a similar shape for z < 0.5, 
since the profile in this region is constrained by the SN 
data, although the underdensity of the spline extends 
to slightly higher rcdshifts. In the polynomial -I- shell 
model, the shell, located at z ^ 2, is quite separate from 
the void, but for the spline model the shell smoothly 
matches onto the void. The presence of positive curva- 
ture is highly significant, with shell amplitude 5i > at 
8cr for polynomial -|- shell. Similarly, the polynomial -t- 
curvature model has ^i > at 7a, and density parameter 
n\°'^ = 1.146 ± 0.013 evaluated at {to, tlss)- Each model 
fits the 557 SNe quite well, the best-fit Xsnc being 534 
for the spline, 539 for the polynomial + shell, and 541 
for the polynomial -I- curvature. 

A more detailed model comparison between the three 
cases is somewhat hard to quantify, since error bars of the 
SNe are tuned to give a x^ P^r degree of freedom of order 
unity for ACDM. We find a corresponding Xsno = ^31 for 
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FIG. 4: Local density parameter ^^ versus redshift z for 
the spline (top panel), polynomial + shell (middle panel), 
and polynomial + curvature (bottom panel) profiles. The 
grayscale level indicates the relative likelihood in the fit to 
the CMB + SN data. 



ACDM using joint CMB + SN data, with the fit to the 
CMB component similar in the void and ACDM models. 
These results indicate that, although the extra freedom 
of the spline improves the SN fit over the polynomial, 
there is still a slight preference for ACDM. This could be 



somewhat alleviated by introducing more control points 
to the spline, but the improvement in A^^ is not sig- 
nificant for the additional degrees of freedom. Even at 
this stage it is unlikely that a model comparison would 
significantly favour ACDM over a three-point spline. 

There are, fortunately, several other parameters which 
are good discriminators between void models and ACDM. 
We show a selection of marginalized likelihoods for the 
polynomial -I- shell profile in Fig. 5, with similar results 
for the spline and polynomial -|- cuvature (the exception 
for curvature being the Compton y-distortion, which we 
discuss later). A general feature of void models is that 
they require an extremely low local Hubble rate, as first 
noted in ZMS. We find /iq = 0.45 ± 0.02 for both the 
spline and polynomial -I- shell, and /iq = 0.47 ± 0.02 for 
polynomial + curvature. As mentioned in ZMS, local 
measurements of h^ (i.e. those independent of the cos- 
mological model) are generally much higher (see, e.g., 
Rcf. [50] for a recent review). For example, recent mea- 
surements of Cepheids in SN hosts at z < 0.1 give 
ho = 0.738±0.024 [51], which rules out the void models at 
high significance. These distances are small enough that 
the dynamics of the void would not significantly affect 
the local distance ladder. Even the lower local determi- 
nations of ho in Ref. [52], who found Hq = 0.623 ± 0.05, 
are strongly at odds with the predictions of the void mod- 
els. (See [53] for a discussion of the reasons for this lower 
ho measurement.) 

The low ho also contributes to what we refer to as the 
"old age" problem for voids. In Fig. 6 we plot the look- 
back time for samples from the MCMC chains. These 
models are significantly older than ACDM, whose age is 
13.7 ± 0.1 Gyr [34]. We find ages 18.8 ± 0.5 Gyr for the 
spline, 18.6 ± 0.5 Gyr for the polynomial + shell, and 
17.6 ±0.4 Gyr for the polynomial -|- curvature. The look- 
back time of void models is already equal to the age of 
ACDM by z « 1. This could pose another serious prob- 
lem for voids, since the existence of observed structures at 
redshifts z > 5, which corresponds to look- back times of 
17-18 Gyr in void models, could be hard to reconcile with 
the ages of the oldest known objects in the Universe, i.e. 
globular clusters, which are consistent with ACDM [55]. 
Although it is difficult to provide an upper limit to the age 
of the Universe from local observations, it is clear that 
there is a serious problem with void models — in such a 
model we should be observing the epoch of formation of 
the oldest known stars at z « 1! 



2. Big bang nucleosynthesis 

Observations of elemental abundances provide a con- 
straint on the baryon-photon number density ratio, j] = 
iT-b/n-y, via our theoretical understanding of the epoch of 
big bang nucleosynthesis (BBN) [56]. For FLRW models 
77 is constant in time, but this is not generally true for 
LTB. (We assume, however, that 77 is spatially constant 
at early times.) Therefore we should explicitly calculate 
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FIG. 5: Selection of marginalized likelihoods for the polynomial + shell profile using the CMB + SN data. For the Compton 
y-distortion the dashed curve shows the addition of radial BAO data, while the vertical line is the CORE 2(7 upper limit [54]. 



r] for our void model near the LSS (or early enough that 
curvature is unimportant). Since the effective EdS model 
has the same early physics at the LSS as our void model, 
we can readily show that the early-time value for voids 
is 
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(41) 



Recent estimates of the deuterium abundance in metal- 
poor damped Lyman-a systems, for example, imply that 
7710 = 5.8 ± 0.3 [57]. From our MCMC chains we obtain 
7710 = 6.1 ± 0.2 for the spline and 6.0 ± 0.2 for the poly- 
nomial -|- shell and polynomial -I- curvature. These are 
entirely consistent with the BBN constraint, which is not 
surprising since the physics at the LSS in the void model 
is the same as that of standard ACDM, and observations 
of the CMB imply a similarly consistent value of 77 within 
the standard ACDM framework. 



3. Dipole 

In Fig. 7 wc show the dipole /3(z) generated from 
our MCMC samples. The peak values are significantly 
higher than our local motion with respect to the CMB, 
for which (3 ^ 10~^ [58]. For small distances from 
the centre of the void, one can write /3 as a linear 
function of the proper distance d. For the spline and 



polynomial (shell and curvature) models we find /3 = 
(3.8 ±0.4) X 10-5dMpc"\ (3.6 ±0.4) x IQ-^dM^c"^, 
and (3.7 ±0.4) x 10-^dMpc"\ respectively. This leads 
to a nominal constraint on how close we must be to the 
centre of the void of ~30 Mpc. However, a more detailed 
analysis including the stochastic (i.e. peculiar velocity) 
and LTB background components leads to weaker con- 
straints [11]. 

The dipole plots show several interesting features, 
which are more transparent in the polynomial + shell 
case. For z < 2 we have Zout > -Zin, so the outgoing pho- 
tons (those passing though the void centre) receive more 
rcdshift than the incoming ones. The initial peak at low 
z results from the higher radial expansion rate inside the 
central void — outgoing photons spend more time in the 
void so are rcdshifted more. There is then an extended 
plateau to z '--^ 2, even though Vt^°^ ^ 1 in this region 
(see Fig. 4). This plateau is due to the Sach- Wolfe (SW) 
effect — outgoing photons originating from the LSS are in 
the region of the overdense shell (on the opposite side of 
the void) , and so redshift out of the gravitational poten- 
tial. Moving further out in redshift, at z > 2 wc have 
^in > ^out- This region results from the lower expan- 
sion rate in the shell, with outgoing photons spending 
more time in the shell. There is also a SW effect from 
the central void around z ~ 3, but this is subdominant 
compared to the blueshifting induced by the shell. 

The general features arc the same for the spline profile, 
but the extended SW plateau cannot be resolved. Im- 
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FIG. 6: Look-back time versus z for the spline (top panel) and 
polynomial + shell (bottom panel) profiles. The polynomial 
+ curvature case is very similar. The dotted curve is the 
look-back time for the WMAP7 ACDM cosmology. 



portantly, for the polynomial + curvature model there is 
no dipolc induced by the homogenous background cur- 
vature. Here the main feature is the peak at low z, with 
the smaller SW signal from the central void also visible 
around z ^ 3. 
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^. Compton y-distortion 

The calculation of the y-distortion in void models 
is complicated by several details. First, the integrand 
in (22) is dependent on j3{z) and hence (for the case of 
zero background curvature) the details of the shell. The 
position of the shell is not constrained by CMB data, 
since the fit to the CMB is only dependent on the inte- 
grated shell density. Second, the integration limit in (22) 
is the reionization redshift Zro, which is uncertain in our 
MCMC analysis. Finally, calculation of y using only the 
dipole calculated from radial geodesies becomes less ac- 
curate in the peripheral region of the void, leading to an 
error on the value of y. 



FIG. 7: Dipole /3 versus z for the spline (top panel), poly- 
nomial H- shell (middle panel), and polynomial + curvature 
(bottom panel) profiles. 



We address these issues by first considering the best-fit 
polynomial -I- shell model. This has y = 3.0x 10~^, which 
is independent (to the first decimal place) of the choice 
of Zro in the range 8 to 15. The best bound on the y- 
distortion comes from the FIRAS instrument on COBE, 
which gives y < 1.5 x 10~^ at 2cr confidence [54]. It is 
clear that this limit could provide important constraints 
on void models. However, on removal of the outer shell 
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(by manually setting 6i = 0) we find i/ = 0.1 x 10~^, 
again independent of Zioi which is a significantly weaker 
constraint. This then suggests that by pushing the shell 
sufficiently far outwards (i.e. past Zrc), we may be able 
to reduce y to satisfy the COBE constraint while still 
fitting the CMB (by adjusting the shell's amplitude — 
recall, however, that we must maintain K{r) < 1 and 
avoid shell crossings). Alternatively, one can evade the 
COBE constraint by considering models with non-zero 
background curvature. Here, the integrated I3{z) is much 
smaller, as shown in Fig. 7, with a best fit of y = 0.2 x 
10~^. From our MCMC chains in this case, we find y < 
2.2 X 10~^ at 2fT confidence, i.e. the majority of samples 
are within the COBE limit. 

We investigate the polynomial + shell model further 
by computing the y-distortion for a grid of shell profiles 
(constrained to have the Tq'^^ and Hq*^^ values required 
to fit the CMB) for the best-fit polynomial central void, 
but we were unable to find a profile which substantially 
reduced y. Shells close to the LSS must have a small 
comoving thickness, and these violate K{r) > 1. The 
most distant shell which still fits the CMB is located at 
z ^ 20, but the dipole for this model has an extended 
SW plateau to z ~ 10 and the resulting y-distortion is 
similar to the nominal best fit. The minimal y is found 
for a shell located at z '--^ 5 with a thickness of 6z ~ 4. 
This has y = 1.0 x lO"'^ with z^o = 8, and y = 1.5 x 10"^ 
with Zj-c = 15. A detailed calculation for this model, 
using non-radial geodesies [37], finds the actual value of 
y to be around 10% higher than the dipole-approximated 
value. 

With these uncertainties in mind, the y values from the 
MCMC chains of models with shells, which we show in 
Fig. 5, could be a factor of 2-3 lower if the shell position 
were carefully tuned. The COBE result then puts mod- 
els in the tail of the y distribution under severe tension. 
These models correspond to wider, deeper voids with a 
larger dipole term, which can be seen in Figs. 4 and 7. 

Finally, note that in Ref. [2] it was found that very 
large voids without overdense shells were already in con- 
fiict with the COBE y-distortion constraint. Thus it may 
seem surprising that our models with large shells, and 
their consequent dipole contribution (cf. Fig. 7), only 
marginally exceed the COBE limit. Part of the expla- 
nation appears to be that Ref. [2] imposed ho ~ 0.73, 
while we have found that much lower values of ho are 
required to fit the CMB. Our lower ho implies a lower 
density of scatterers and hence lower y-distortion. In 
addition, Ref. [2] employed a simplified linear "Hubble 
bubble" model, rather than performing a rigorous LTB 
treatment. 



5. Radial BAO 

Estimates of the radial BAO (RBAO) scale have been 
made recently in Refs. [59, 60]. These were used in the 
ZMS analysis (see also [22]), but since that time the sys- 
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FIG. 8: Radial BAO scale Az versus z for the spline (top 
panel) and polynomial -I- shell (bottom panel) profiles. The 
polynomial + curvature case is very similar to the polynomial 
+ shell, but without the suppression at z ~ 2. Also shown 
(dotted curves) is the best fit ACDM model, along with the 
measurements from Refs. [59, 60]. 



tematic errors have been revised and increased by a factor 
of about a third. In addition, the statistical significance 
of the claimed RBAO detection has been questioned in 
Refs. [61, 62] (however, see [63] for a counter-argument). 
With this caveat in mind, in Fig. 8 we show the RBAO 
scale Az for our MCMC chains, along with the best fit 
ACDM model and the measured Az at z = 0.24 and 0.43. 

One can immediately see that Az is extremely sensi- 
tive to the particular profile. The polynomial -I- shell 
models have two regions of strong suppression of Az — 
one on the edge of the central void, where matter piles up 
after fiowing outwards from the void (this can be seen by 
closely inspecting Fig. 4), the other corresponding to the 
outer shell. In the spline models these two regions are 
indistinguishable. This suppression of Az is due to the 
low radial expansion rate in those regions, as explained in 
Sec. Ill A 5. The polynomial -I- curvature model is very 
similar to the polynomial -I- shell, without the suppres- 
sion at z ~ 2. 

Post processing our MCMC chains with the addi- 
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tional RBAO data increases the x^ of the best fit by 
Ax^ = 10.1, 7.8, and 7.5 for the sphne, polynomial + 
shell, and polynomial + curvature cases, respectively. 
The large increase in x^ shows that current RBAO data 
are already a strong discriminator from ACDM. More- 
over, RBAO prefers wider void profiles, in opposition to 
the y constraint. In Fig. 5 we show the marginalized y 
values for the polynomial + shell model with the addition 
of RBAO data. Even with factors of 2-3 uncertainties in 
y (from tuning the position of the shell), voids with shells 
are strongly disfavoured considering the COBE limit. 

For the polynomial + curvature model we find the 
posterior is approximately Gaussian, with y = (1.4 ± 
0.7) X 10"'^, when including RBAO. Although closer to 
the COBE limit, the wider voids favored by RBAO are 
still compatible with the y constraint when curvature is 
distributed homogeneously, rather than in the form of an 
overdense shell. 



6. Local (jg 

The marginalized values of the matter power amplitude 
CT8 near the void centre today are 0.48 ±0.04, 0.49 ±0.04, 
and 0.48 ± 0.04 for the spline, polynomial + shell, and 
polynomial + curvature, respectively (see also Fig. 5). 
These values are much lower than that in ACDM [34]. 
The main reason for this is that near the centre the LTB 
model is well approximated by an open FLRW spacetime, 
as explained in Sec. Ill A 6. For fixed Jim today, the 
growth of structure is much more strongly suppressed in 
open as opposed to A models (see, e.g., [64]), leading to 
lower cTg and interesting constraints. As with the low Hq 
and old age problems already discussed, such low local ag 
values appear to be a generic problem with void models. 

Local measurements of ag include those from weak 
lensing, the galaxy power spectrum, and cluster abun- 
dances. In order to compare our results with observa- 
tions, we should ensure that the data analysis was not 
model-dependent, for example by implicitly assuming a 
ACDM background (as the CMB-based measurements 
of as do). One such constraint comes from the mass 
function of local galaxy clusters [65]. There is a strong 
correlation between the cluster number density and erg, 
and simulations show a similar functional form in both 
open and ACDM models [66]. The complication comes 
from relating the cluster mass with an observable prop- 
erty such as temperature or luminosity. 

Recently, there has been significant progress in un- 
derstanding systemic errors which affected previous con- 
straints. Ref. [67], for example, provide a detailed dis- 
cussion of the various sources of error. They use a local 
sample of clusters (46 of the 48 objects are at z < 0.1, 
the other two at z < 0.2), and find as (fiin/0.32) ' = 
0.86 ± 0.04 (for fin, < 0.32) by relating the cluster count 
to X-ray temperature. The error budget is dominated by 
uncertainties in cluster physics, and cosmology depen- 
dence is negligible at such small z. A summary of other 



recent as results can be found in Ref. [68]. 

For void models, with the Jlni = ^m'^(^ = 0) likeli- 
hood peaked at « 0.2 (cf. Fig. 5), the low value of as we 
have found is incompatible with this limit, at very high 
significance. 



IV. EXTENSIONS 

A. Modifications to the initial power spectrum 

For a simple power-law primordial spectrum and a spa- 
tially fiat background (as the simplest models of infia- 
tion predict), voids without an outer shell do not fit the 
CMB at many a, since the physics at the LSS is so dif- 
ferent from ACDM. We found that this result persists 
even when we allowed for the running of the spectral in- 
dex, Eq. (34). Intuitively, the redshift from the void is 
small compared to that from the shell, so without a shell 
we expect T^'^^ « Tq, whereas recall from Sec. IIIAl 
that we need T^'^^ sa 3.4 K to fit the CMB. One way to 
try to compensate for the different physics is to consider 
non-trivial modifications of the primordial spectrum. A 
model which has been discussed in the literature (see 
e.g. [69]) is a broken power law (BPL). This was origi- 
nally proposed to try to make EdS compatible with ob- 
servations without the need for dark energy. The four- 
parameter spectrum has the form 



Ps{k) 



As{k/kor-\ 



k < ki 



k> ki, 



(42) 



where ko is the normalization scale and a matching con- 
dition in Vs{k) is enforced at the break scale fci. The 
spectral index is ni for k < ki and n2 for k > ki. 

To test whether we could obtain an improved fit with 
this spectrum, we repeated the effective EdS analysis 
with the BPL spectrum. We fix T^'^^ = Tq since this 
is a feature of voids without shells. The best-fit likeli- 
hood relative to ACDM is Ax^ = —14.2, and we show 
the CMB spectrum for each model in Fig. 9. The agree- 
ment at low i is good, but the lack of power at high £ 
results in a very poor overall fit with a large Ax^ (note 
that the use of the small-scale CMB data was impor- 
tant here). In both cases we do not include lensing of 
the CMB, but since lensing only smooths out the acous- 
tic peaks we do not expect this to alter our conclusions. 
This means that current CMB data are already sensitive 
enough to disfavour EdS with a BPL and T^'^^ = Tq 
over ACDM. Hence further features in the primordial 
spectrum would need to be introduced in order to try 
to improve the fit. The BPL marginalized parameters 
arclogio[fci-Mpc] =-1.90 ±0.04, with rii =0.77 ±0.01 
and 7i2 = 1.33±0.08. It is worth stressing that this corre- 
sponds to a very strongly but oppositely tilted spectrum 
on either side of the break, i.e. to a very large depar- 
ture from the spectra that the simplest models of infia- 
tion predict. The implication for voids is that shells or 
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FIG. 9: High (. CMB spectrum. The best-fit ACDM model 
is shown by the sohd curve, the effective EdS with power- 
law spectrum by the dotted curve, and the effective EdS with 
fixed T^ — To and a broken power law by the dashed curve. 
Also shown is binned WMAP and ACBAR data. 



background spatial curvature, and their consequent very 
low Hq, are still required in order to fit the CMB, even 
with a four-parameter BPL. Furthermore, any attempt 
to improve the fit without shells or curvature would en- 
tail considerable tuning of the primordial spectrum and 
large departures from scale invariance. 



B. Multi-valued profiles 

1. Parameter space 

Multi-valued solutions arise when the radial expansion 
rate Y' /Y' , and hence, by Eq. (9), dz/dr, becomes neg- 
ative. This can occur in overdense regions, either from 
matter piling up after flowing outwards from the void, 
or within the outer shell. In this regime there are three 
solutions to r{z), and hence three "branches" of dp^ or jjl 
for a given redshift. 

For the polynomial -I- shell profile the multi-valued re- 
gion of parameter space is far from that found using the 
CMB -I- SN constraints in Sec. IIIC. The spline allows 
us to investigate wider voids with a sharper transition 
from the void to overdense shell. These models have a 
multi-valued region closer to that allowed by the data, so 
we focus on these in the following analysis. 

We first perform a survey of the parameter space. To 
do this we modify our code to use {rg = 0,ri,r2,52 = 
0,f]JS'=(z = 0),ro^'*s,/i|^'iS} as the input parameters 
rather than {tq = 0, ri,r2,(5o, 5i,(52 = 0, /iq} via an it- 
erative search method. We then fix T^'^'^ and h^"^^ to 
their best-fit values from Sec. IIIBl, so the profiles fit 
the CMB. In Fig. 10 we show the reduced x^ from a fit 
to the Union2 SN data for various values of n|°'^(z = 0). 
We also indicate the value of 100 h^ and the region where 
multi- valued solutions exist. For multi- valued profiles we 



compute the x^ in the "kindest" possible way — that is, 
we compare the observed value with the closest branch 

of /i. 

For single-valued solutions it is clear that deep voids 
[i.e. models with low Vl^°^^{z = 0)] do not improve the fit to 
the data. However, these deep models have an interesting 
multi-valued region, where the "kind" y^ is significantly 
better than the single-valued solutions. An example of 
one of these profiles is shown in Fig. 11. The low y^ value 
results from improved fitting of SN outliers around the 
branches with higher ^. However, most of the samples 
are closest to one branch, which has the lowest /i. If this 
was indeed the true underlying model, we might expect a 
larger fraction of SNe around the upper branches. An im- 
proved x^ calculation should therefore take into account 
the expected distribution of SNe around the branches. 



2. Modified supernova fit 

Our modified x^ consists of two components: (1) Xkind' 
which we mentioned above; and (2) Xdisf from the ex- 
pected distribution of SNe around each branch. For Xdist ' 
suppose there are TV supernovae in the multi-valued re- 
gion at redshifts Zi (i = 1,...,N) with three magni- 
tude residuals fii{zi) < ^2{zi) < ^^{zi), as illustrated 
in Fig. 12. Each magnitude residual /x^ has standard de- 
viation 5i. 

Now, suppose all supernovae are distributed around 
branch 2 according to a Gaussian probability function. 
Then the probability of the i'th supernova lying closest 
to branch 1 or 3, but being associated with branch 2, is, 
respectively. 



P21, 
P23, 



1 — erf 



1 — erf 



/^2 - A^i 

fJ-3 - tJ-2 

2V2S^ 



(43) 
(44) 



where erf is the error function. For branch 2, the proba- 
bility is P22,i ~ 1 —P2i,i — P23,i- The cxpected number of 
supernovae Nj closest to each branch j, along with the 
variance cr^ , can be derived by considering multiple trials 
from a Bernoulli distribution. The result is 



A^, 



N 

E 



N 



P2j, 



'^P2j.^il~P2j,^) ■ (45) 



Next, we allow branch 1 or 3 to be the "correct" model. 
Similar expressions can be derived for the probability of 
each supernova lying closest to a branch, but being asso- 
ciated with another. Finally, we suppose that supernova 
i has a probability of belonging to each branch j. We 
label this branch probability bj^i, where '^j^ibj.i = 1. 
We assume that the branch probability is proportional 
to the number density of SNe per redshift interval dz at 
that redshift, such that 



ex 



dN _ pY^Y' 
dzdn ~ {l + z)Y'' 



(46) 
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FIG. 10: Multi- valued region of parameter space for the spline profile. The top left panel has r2m'^(2 = 0) = 0.2, top right 0.15, 
bottom left 0.1, and bottom right 0.05. Shaded contours show the reduced x from a fit to SN data (note the different scale 
in the bottom right panel) , and dashed contours the value of 100 ho . The region at the bottom of each panel below the heavy 
solid line is the multi- valued regime. 



In doing this we are implicitly assuming that the survey 
is equally sensitive to each of the branches, and that there 
is no SN evolution. The expected number of SNe closest 
to each branch is then 



AT 3 

^J =5m^'='*^'=J'■'• 
^=l fc=l 



(47) 



The variance calculation is more complicated, but we can 
use the approximation 



N 



'^^ ~J2J2p''^'''^^^p''^''^ ■ 



(48) 



-t=l k=l 



This quantity is larger by a factor of approximately 1- 
2 than the actual variance (we have checked this using 
numerical realizations), so is a conservative estimate of 



cr|. The "distribution" component of the x^ is then 

3 /,x ,r„K=\2 



Adist 



E 



[N, - Nf^ 



(49) 



where N°^^ is the actual number of SNc closest to each 
branch. 



After modifying the x^ according to this procedure, 
multi- valued solutions no longer give an improved fit. In 
other words, the "kindest" approach to x^ was too kind. 
We investigated the multi-valued parameter space for 
deep voids over the range rjj^'^(z = 0) = 0.02-1 and found 
that the best fit had a similar x^ jN w 1 to the single- 
valued models. A more detailed study of the branch prob- 
abilities and magnitude errors in these models could be 
interesting and possibly offer some improvement over the 
fit to the SNe of "regular" voids. 

However, the deep voids still suffer from the same in- 
herent problems in fitting other cosmological data which 
the single-valued models exhibit. In particular, they do 
not fit the CMB without an extremely low Hq, as Fig. 10 
shows. Note that this contradicts the results of [70], who 
found that extremely deep voids could fit the CMB + 
SN data with reasonable Hq , although they stressed that 
such models were implausible. The multi- valued regions 
in deep models are also very close to shell-crossing sin- 
gularities. Therefore it is important to treat LTB mod- 
els with multi-valued regions with great care. Finally, 
we note that a dense spherical shell of matter at z < 1 
would likely leave other observable signatures on struc- 
ture which we have not considered here. 
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FIG. 11: Example of a multi-valued profile for ilj^'^(2; = 0) = 
0.05. The top panel shows the density parameter and the 
bottom panel shows the SN magnitude residual, together with 
the Union2 data. The "kind" x^ is significantly better than 
that of the single- valued solutions. 



V. CONCLUSIONS 

We have performed the most complete and up-to-date 
analysis of the proposal to explain the acceleration of the 
Universe with a large, nonlinear void. Our main results 
can be summarized as follows: 

• Models which fit both the CMB and SNe have an 
extremely low local Hubble rate of ho sa 0.45 ± 0.02 
(with the precise constraint being mildly model de- 
pendent) and are very old, with an age of around 
17.5-19 Gyr. Both of these are inconsistent with 
observations. The use of the full CMB spectra, 
rather than just an acoustic angular scale, was es- 
sential for obtaining this result. 

• The radial BAO scale is strongly suppressed in the 
outer void region compared to ACDM, and the data 
are a poor fit to the prediction for voids. 

• The Compton y-distortion constraint from COBE 
rules out wider, deeper voids (with shells) even if 
the overdense shell is tuned to minimize y. Generic 
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FIG. 12: Illustration of SNe distributed in the multi-valued 
region of the redshift-luminosity distance relation. The hori- 
zontal axis is a measure of redshift. 



shells are inconsistent with the COBE limit and 
all void models (with shells) are ruled out at high 
significance in conjunction with radial BAO data. 
Models with homogenous background curvature 
can, however, evade the COBE constraint. 

• The local amplitude of matter fluctuations is ex- 
tremely low in void models, with as « 0.49 ± 0.04. 
This is due to the suppression of growth compared 
to ACDM, and is strongly incompatible with esti- 
mates from local galaxy clusters. 

• The primordial power spectrum would need to be 
highly tuned and far from scale-invariant in order to 
attempt to alleviate some of these problems. Run- 
ning of the spectral index and a four-parameter 
broken power law were insufficient to change our 
conclusions. 

• We uncovered a class of models with multi-valued 
distance-redshift relations. While unusual and in- 
teresting to study, they again did not alter our con- 
clusions. 

• All of our results, with the exception of the 
y-distortion constraint, persisted regardless of 
whether spatial curvature outside the void was dis- 
tributed in an overdense shell or homogeneously. 

Thus it appears that void models for acceleration are 
overwhelmingly at odds with several types of observa- 
tions. These conclusions are insensitive to the details 
of the void profile chosen, as is apparent from the large 
range of profiles our three parameterizations can capture, 
or to the presence of background spatial curvature. Ei- 
ther a large, overdense shell or substantial positive back- 
ground spatial curvature is needed to fit the CMB, al- 
though both move us even farther from the standard in- 
fiationary scenario than an isolated void in an EdS back- 
ground. 

We have assumed that any LTB decaying mode is ir- 
relevant today, although it is likely that a significant de- 
caying mode could help these models evade some of the 
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constraints we have examined here. However, as dis- 
cussed in Ref. [27], a substantial decaying mode today 
requires the early Universe to be extremely inhomoge- 
neous, which is completely at odds with the standard in- 
flationary scenario. Ideally, though, models with decay- 
ing modes should be confronted with observations rather 
than argued against on such grounds [37]. 

It is worth stressing the crucial role that CMB ob- 
servations have played in our study. The situation has 
a close parallel within the FLRW class of models. As 
is well known (see, e.g., [71]), observations of the CMB 
alone are consistent with closed models with considerable 
curvature (and even vanishing A). However, once local 
observations of the Hubble rate are included, a flat model 
is strongly favoured. Similarly, the CMB alone docs not 
rule out the inhomogeneous-curvature void models, but 
does in conjunction with local measurements of Hq. 

Finally, we note that two of the most important as- 
sumptions in cosmology are those of the cosmological and 
Copernican principles. Therefore, in confronting void 
models, which blatantly violate both of these principles, 
with observations, we do more than just examine an un- 
usual approach to the mystery of acceleration. We put 
the foundations of modern cosmology themselves to the 
test. 



Note Added 

Almost simultaneously with the submission of the first 
version of this paper, two other works appeared address- 
ing similar questions [72, 73]. The conclusions of these 
papers appear to be very difi^erent from ours. For exam- 
ple, the abstract of Biswas, Notari, and Valkenburg [72] 
(hereafter BNV) states "We find that the inclusion of a 
nonzero overall curvature drastically improves the good- 
ness of fit of the Void model, bringing it very close to 
that of a homogeneous universe containing Dark Energy, 
while by varying the profile one can increase the value of 
the local Hubble parameter which has been a challenge for 
these models''. 

We agree with BNV that a near-compensated void (i.e. 
a void without an overdense shell) in an EdS background 
cannot fit the CMB regardless of the value of Hq, as was 
already clear in ZMS. However, one of the main claims of 
BNV is that the inclusion of non-zero background spatial 
curvature improves the goodness of fit dramatically. For 
this reason we enlarged the parameter space of models 
in our MCMC analysis to include background curvature 
(instead of overdense shells). Our results, as reported 
above, are essentially unchanged from the overdense shell 
case. In particular, it does not matter how the "positive 
spatial curvature" required to fit the CMB is distributed, 
whether in the form of an overdense shell or homogeneous 
background curvature: the local Hubble rate Hq is too 
low in void models. Importantly, we agree with BNV that 
the presence of background curvature (or an overdense 
shell as we found in ZMS) does dramatically improve the 



fit to the primary CMB anisotropics. However, we still 
find Hq values unrealistically low. 

How do the Hq values presented here compare with 
those of BNV? Our results in Sec. HI C 1 correspond to 2a 
upper limits of /iq = 0.50 for the spline parameterization, 
ho ~ 0.49 for polynomial -I- shell, and /lo = 0.52 for 
polynomial + curvature. The best-fit values of h^ for the 
CMB + BAO + SN -h HST constraints in Tabic 13 of 
the published version of BNV range from 0.495 to 0.552 
for BNV's profiles A to D. However, we have noticed 
a discrepancy with the values in that table. When the 
background LTB shear is negligible, the LTB evolution 
must match that of a spatially curved FLRW model, as 
we mentioned in Sec. Ill A 6. This means that at the 
origin, as well as asymptotically outside the void, the 
LTB models must satisfy the FLRW consistency relation 



Ht = 
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''^^, cosh-i| 
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K^ 


\\-Vt\- 



(50) 



where ^^ = 1 — ^^ , and each quantity is evaluated 
at the same spacetime event. We have noticed that the 
"in" values for profiles A to D in the published version 
of BNV violate the relation (50) at the 5-10% level. If 
the BNV Hq values are adjusted to satisfy Eq. (50), they 
must be lowered by 5-10%, which would bring them in 
line with our 2a upper limits. Of course, it is not clear 
in the first place that the 2a upper limit of a fit to CMB 
-I- SNe (as we have done) is expected to agree well with 
the best fit to CMB -f- BAO + SN -(- HST (as BNV have 
done). 

In addition, we point out that there exists a relation 
between the central density parameter and Hubble rate, 
such that somewhat larger Hq values can be obtained 
with larger fi'^'^g values, as our Fig. 10 shows. BNV's 
profile D obtains a slightly higher Hq by pushing J^'^'^q 
very high (the original version of BNV contained profiles 
with extremely large f^'^'^o)- Nevertheless, we find that 
even by pushing i^'^'^g to unrealistically high values [48, 
49], we cannot significantly increase our upper limits on 
Ho (see Fig. 10). 

We have also attempted to recreate BNV's profile E. 
To calculate the CMB for this model, we must create a 
(closed) FLRW model that shares the same local den- 
sity parameter at redshift z^ as the (closed) LTB model. 
However, we find that the proper circumference of the 
LTB model at z^ exceeds the equatorial circumference of 
the FLRW model by tens of percent, so that the matching 
is not possible. What this tells us is that profile E does 
not approach homogeneity sufficiently rapidly outside the 
void in order to reliably generate an effective model to 
calculate the CMB. Hence we cannot quantitatively trust 
the results for profile E. 

Finally, we note that BNV claimed that their mod- 
els with the largest Ho were compatible with the local 
Hubble parameter measurements in Ref. [74], namely 
ho = 0.623 ± 0.063; which are both lower and have much 
larger error bars than the more recent result from [51], 
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who find ho = 0.738 ± 0.024. Therefore, as BNV point 
out themselves, the models they examined exhibit sub- 
stantially worse fits to the newer Hubble rate data. 

To summarize, we find the largest ho in void models 
that provide a good fit to CMB + SNe to be around 
0.49-0.52. It appears that the results of BNV are es- 
sentially consistent with ours, so that the Hubble values 
in void models arc strongly at odds with the latest local 
measurements . 

The second paper that appeared almost simultaneously 
with ours was Ref. [73], whose abstract states "We allow 
for the dynamical effects of radiation while analyzing the 
problem, in contrast to other work which inadvertently 
fine tunes its spatial profile. This is a surprisingly im- 
portant effect and we reach substantially different con- 
clusions. " Essentially, the authors of Ref. [73] claimed 
that allowing for an 0(1) radiation- matter isocurvature 
mode at early times can provide a loophole to our re- 
sult that the CMB cannot be fit without a very low Hq. 
First of all, we point out that the matching procedure we 
use to generate the effective EdS model does not ignore 
radiation at early times; rather, it ignores the shear for 
2 > Zni in effectively matching to a FLRW matter + ra- 
diation -|- curvature model at early times. Our method 
uses the LTB solution (with test radiation source) at late 
times when that is a very good approximation, and uses 
FLRW at early times when that is a good approximation. 

We maintain that an 0(1) radiation-matter isocurva- 
ture mode at early times cannot have an important ef- 
fect at late times. To see this, recall that observations of 



the CMB by the central observer specify the physics at 
the LSS, i.e. the matter and radiation densities, via the 
shape of the Ci anisotropy spectrum. Therefore the pres- 
ence of any such isocurvature mode would imply matter- 
radiation ratios different from the standard ACDM values 
at early times inside our past light cone. But any such 
excess or deficit of radiation would free stream away and 
could not affect the central observers today, since by con- 
struction it does not lie on their past light cone. At worst, 
the late time central effect of an early central radiation 
over- or underdensity can only be a decaying LTB mode, 
corresponding at late times to an effectively inhomoge- 
neous bang time, tB{r). But the amplitude of such a de- 
caying mode could be at greatest O(ioq), where tcq is the 
time of matter-radiation equality. Therefore, the effect of 
such a decaying mode today, which would be determined 
by the ratio t^q/to, would be negligible. Nevertheless, we 
believe that a more thorough examination of the effects 
of radiation in void models is warranted [37]. 
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